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ABSTRACT 

We investigate the effect of star formation and diffuse photoelectric heating on the properties of 
giant molecular clouds (GMCs) formed in high resolution (< 10 pc) global (~ 20kpc) simulations of 
isolated Milky Way-type galaxy disks. The clouds are formed through gravitational fragmentation 
and structures with densities njj, c > 100 cm~ 3 are identified as GMCs. Between 1000-1500 clouds 
are created in the simulations with masses M > 10 5 M Q and 180-240 with masses M > 10 6 M Q in 
agreement with estimates of the Milky Way's population. We find that the effect of photoelectric 
heating is to suppress the fragmentation of the ISM, resulting in a filamentary structure in the warm 
gas surrounding clouds. This environment suppresses the formation of a retrograde rotating cloud 
population, with 88% of the clouds rotating prograde with respect to the galaxy after 300 Myr. The 
diffuse heating also reduces the initial star formation rate, slowing the conversation of gas into stars. 
We therefore conclude that the interstellar environment plays an important role in the GMCs evo- 
lution. Our clouds live between — 20 Myr with a high infant mortality (t' < 3 Myr) due to cloud 
mergers and star formation. Other properties, including distributions of mass, size and surface density 
agree well with observations. Collisions between our clouds are common, occurring at a rate of ~ 1/4 
of the orbital period. It is not clear whether such collisions trigger or suppress star formation at our 
current resolution. Our star formation rate is a factor of 10 higher than observations in local galaxies. 
This is likely due to the absence of localized feedback in our models. 

Subject headings: galaxies: spiral, galaxies: ISM, galaxies: star clusters, methods: numerical, ISM: 
structure, ISM: clouds, stars: formation 
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1. INTRODUCTION 

Embedded in the turbulent gas of the interstellar 
medium (ISM) are the stellar nurseries of the Galaxy; 
the giant molecular clouds (GMCs). The properties of 
these clouds dictate the environment in which most new 
stars are formed, making these objects the primary con- 
trollers of star formation in the Galaxy. To understand 
star formation, it is therefore essential that we first un- 
derstand the processes that govern the formation and 
evolution of the GMCs. 

Exactly what determines the GMC properties is a 
taxing problem. Theoretical work on their forma- 
tion has led to two different schools of thought. The 
"Top-down" mechanism suggests GMC formation is 
driven via large-s cale gravitational or magnetic disk 
instab ilities (e.g. IShettv fc Ostrikerl 120061: Kim et al.l 
120031 iGlover fe Mac Low! l2007al |b[). whereas "bottom- 
up" processes see c louds formed in colliding flows (e.g. 
IHeitsch et all [20081 or via agglomeration from inelastic 
collisions between the GMCs (|KwanllT979D . It is possi- 
ble that both these pr ocesses are im portant in different 
galactic environments ([Dobbs 20081). 

Depending on the lifetime of the cloud, the impor- 
tance of the formation mechanism on the GMC prop- 
erties varies. If the cloud lives as long as a few free-fall 
times, then its attributes could depend primarily on in- 
teractions with its environment (e.g. disk sheer, spiral 
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arms and cloud-cloud collisions) , and also on its internal 
processes driven by star formation, including turbulence 
from supernovae and radiation pressure. Whether GMCs 
are long-lived enough for this is highly debated with the- 
ories supporting lifetimes both short and long compared 
to the free-fall timescale, 
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However, current estimates s uggest that a GMC lives 
between 1-2 free-fall times (e.g. iMcKee fc Ostriker|[2007i 
and references therein). This implies that the cloud's 
evolution could play a significant role in determining its 
properties. 

Observationally, it is difficult to measure the cloud 
properties accurately since the molecular gas cannot be 
measured direct ly, but must be inf e rred from the abun- 
dances of CO (IGlover et al.1 [2010t IShettv et al.ll20Toh . 
From the theoretical stand point, there is the prob- 
lem of scale that the forces are operating on. While 
the GMC itself is of order 20 pc in radius, it lives in 
a galactic disk whose stellar component, in the Milky 
Way, extends to around 20kpc in radius. This means 
to self-consistently replicate their properties, simulations 
of GMCs must encompass 3-4 orders of magnitude in 
scale. Previous work that has studied the properties of 
GMCs on the g alactic scale has often b een limited to 
two dimensions ()Shettv fc Ostrikerl 12008) or has had to 
assum e a fixed two-phase medium for the ISM (|Dobbsl 
2008). While the majority of the clouds are confined to 
the plane of the disk, both cloud collisions and feedback 
eject gas from the surface and, in the latter instance, 
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it is this vertical expulsio n that is thought to regulate 
the pressure in the ISM (|McKee fc Ostrikeii H977t ICoxl 
2005); the environment in which the clouds are forming. 
Similarly, previous three dimensional studies of the ISM 
in galaxy disks at lower resolution, produce a continu- 
ous multiphase medium tha t is poorly represented by a 
fixed, discrete phase m odel (jTasker fc Bry an 2006, 120081: 
iWada fe Nor^nar][200l . 

Local models which consider a partial region of the 
galaxy disk (normally less than a few kpc across) are 
able to achieve much higher resolutions and -with the 
use of a shearing box- can approxim ate the effects of 
galactic shear on the cloud proper ties (|Kim et alJ 120081 : 
iKim fe Ostriken 120071 120061 l2001ft . Alternative models 
on si milar scales investigate th e results from colliding 
flows (|Heitsch et al . 2008, 2009) and the local impact of 
supernovae and t urbulence fe.g. Uoung fc Mac Lowj|2006l ; 
ISlvz et aT]|2005l ). These simulations are able to tell us 
much about the structure of the GMCs and star-forming 
gas, but they cannot explore the evolution of the GMCs 
as they move through the disk or compare with globally 
averaged properties. 

Whether the global environment is important for cor- 
rectly modeling the evolution of the GMCs remains an 
open question. If GMC properties are not a strong func- 
tion of their environment, then they can be modeled as 
separate entities and the inclusion of global forces from 
the galactic disk are not necessary. This hypothesis is 
given some support by the observational results from 
populations of GMCs in galaxies other than the Milky 
Way, which find their properties are similar to those of 
Galactic GMCs, including their velo city dispersions, sur- 
face density and virial parameters ( [Bolatto et alJ [20081: 
iRosolowskv et alJ 120031 : lFukui et al.ll2008l ). IBigiel et alJ 
(2008) also found that the 18 galaxies in The HI Nearby 
Galaxy Survey (THINGS) appeared to have a fixed star 
formation rate (SFR) per unit of molecular gas, all of 
which could suggest a set of universal properties for 
GMCs. 

On the other hand, if the clouds are gravitationally 
bound, then their velocity d ispersions become e stab- 
lished via GMC interactions (|Gammie et~aTlfl99l . In 
this case, their radial position in the galaxy can change, 
causing greate r susce ptibil ity t o gala ctic shear. Both 
iGammie et ail (|1991ft and [jai3 (|2000l) argue that self- 
gravitating GMCs should suffer relatively frequent col- 
lisions, which could be an important process in control- 
ling the molecular mass of a GMC. This would impact 
the GMCs SFR, making the galactic environment an es- 
sential ingredient in the GMCs evolution. Furthermore, 
while the mass profile of GMCs is universally seen to be 
a power law of the form, 



the value of t he exponent, a c , is f ound to differ be- 
tween galaxies. IWilliams fc McKeel (|1997h measure a 
value of a r between 0.6 — 0.8 in the Milky Way, while 
Rosolows kv et all (I2003D finds a st e eper g radient of a c ~ 
1.6 in M33. iBlitz fc Rosolowskv! (|2004h conclude that 
these variations are not due to systematic uncertainties, 
but are a consequence of the galactic environment. This 
suggests that the global-scale structures impact the range 



of masses of clouds formed. 

Perhaps most suggestive evidence in favor of an in- 
trinsic connection between clouds and their g alactic en- 
viron ment is the Kennicutt-Schmidt relation dKennicuttl 
1998). This empirical relation links the averaged gas sur- 
face density, E gas to the surface density of star formation 
rate E s f r via: 

E sfr cx (3) 

where iKennicuttl ()1998l ) found the exponent a s f r = 1.4 ± 
0.15 for spatial averaging over the en tire disk. More re- 
cent studies using the THINGS data (|Bigiel et al.ll2008| ) 
find the relation to hold linearly with the molecular com- 
ponent of the gas with a a s f r = 1.0 ± 0.2, with spatial 
averaging (resolution) of 750 pc. 

An alternatively fit for this data can be found by re- 
lating S s f r to the orbital angular frequency at the outer 
radius, Sl out : S s fr cx ~S g Q OU f This result applies even 
in the more extre me environments of s tarburst and high 
red-shift galaxies (jGenzel et al.l [2010l ). Both these re- 
lations indicate an intimate connection between global 
structure and the GMC star- forming environment. 

On scales of the same order as the GMC size, an im- 
portant component in the GMC evolution must come 
from the stars themselves. Star formation is observed 
to be highly clustered, with star cluster s forming out 
of de nse clumps with initial radii ~ lpc (|Lada fc Ladal 
2003). Within this small region, the total star formation 
efficiency is relatively high at ~ 0.1 — 0.5, but the major- 
ity of the GMC is not for ming stars, poss ibly due to the 
effect of magnetic fields (|Crutcher 1 12005) . The average 
efficiency over the whole cloud is therefore of order a few 
percent per local free-fall time (Krumho lz fc Tanl 120071 : 
IZuckerman fc Evansl[l97l . 

Gas is removed from the cloud to create a star 
which then deposits energy into its surrounding medium 
through diffusive and energetic feedback. During its life- 
time, a massive star will be source of FUV radiation 
which can be absorbed by dust grains to eject an elec- 
tron that will heat the gas. This photoelectric heating 
has long been thought to be the dominant form of heating 
in the neutral ISM, which includes the GMC population 
(jWolfire et all fl995h . More energetic forms of feedback 
from stellar winds and supernovae will also act to deposit 
concentrated blasts of energy into the star's immediate 
environment. Whether t he cloud can survive the star's 
life-cycle is debated (e.g. iMurravlfeOlOl ) but the fact this 
process will affect the cloud's properties is not. 

Due to the range of the forces in play, it is exceptionally 
difficult to determine the dominant processes affecting a 
GMCs evolution. Is it the interaction with the global 
galaxy environment, the results of star formation and 
feedback or an equal combination of systems? In these 
papers, we aim to investigate this with a set of global 
disk simulations that separate out the processes by in- 
troducing each influenci ng factor indiv i dually . 

In our first paper (jTasker fc Tanl 120091 hereafter 
TT09), we simulated an idealized population of GMCs 
without the presence of star formation or any form of 
feedback. Despite the simplicity of the model, we re- 
produced many of the observable properties of measured 
GMC populations including mass surface density, veloc- 
ity dispersion, angular momentum and vertical distribu- 
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tion. In addition to this, we found a typical collision time 
between clouds of ~ 20 % of the local orbital time, in 
agreement with estimates bv lTanl ([2000D . This suggested 
that compressive flows generated in cloud collisions could 
be a dominant mechanism for inducing star formation. 
Such a process has parallels with the local-scale colliding 
flow models which trigger star for mation via compres- 
sive forcing of the turbulent flows (iBaneriee et al.l 120091 : 
iHeitsch et all [20081: iHennebelle et al.l 12008). However, 
for this simulation without star formation and feedback 
processes included in the model, the GMCs could only 
be destroyed through mergers. This meant that it was 
not possible to regulate their evolution, giving a steadily 
more massive population over time. 

In this paper, we extend the model in two new simula- 
tions. The first of these includes star formation without 
feedback and the second contains star formation with 
diffuse feedback from FUV photoelectric heating. We 
compare the properties of the GMC populations formed 
in both models and with the GMCs formed without star 
formation in TT09. Our results will show that the diffuse 
heating reduces the fragmentation of the disk, causing 
clouds to be embedded in a filamentary warm ISM. This 
reduces the initial star formation rate and suppresses the 
formation of a retrograde rotating population of clouds. 
We will see that both our population of clouds continue 
to match many of the observations of GMCs in the Milky 
Way, including the mass profile, size, mass weighted sur- 
face density, vertical distribution and gravitational bind- 
ing. With the inclusion of cloud destruction through star 
formation and mergers, we are able to estimate the av- 
erage lifetimes of the clouds in our populations which 
are found to be largely between — 20 Myr in agreement 
with current estimates. We will show that our SFR is a 
factor of 10 higher than that observed in local galaxies 
and suggest that this is due to the lack of energetic local 
feedback in our models. 

These localized feedback processes, such as supernovae, 
radiation pressure and ionization are expected to play an 
important role in GMC evolution. However, to under- 
stand the determining forces on the GMC population, 
we focus on the two effects of star formation and diffuse 
heating in this paper, leaving additional forces for later 
study. 

Details of our method are outlined in section §2, global 
properties of the disk are presented in §3 and the struc- 
ture of the ISM in §4. §5 will focus on the properties of 
the individual GMCs and §6 will look at the star for- 
mation in the disk. In §7 we will present our conclusions. 

2. NUMERICAL TECHNIQUES 

2.1. The Code 

The simulations presented in this paper were run us- 
ing Enzo; a three-dimensional adaptive mesh refinement 
(AM R) hydrodynamics codeJO'Shea et al.| [200l iBrvanl 
119991: iBrvan fc Normanlll997t) . The AMR technique is 
particularly adept at resolving multiphase fluids such 
as those found in the ISM, due to the natural bound- 
aries between grid cells producing a ccurate resolution o f 
shocks and low numerical mixing (jTasker et al.l l2008h . 
Enzo has previously been used to model galactic disks 
where it successfully produced a self-consistent atomic 



multiphase ISM dTonnesen fc Bryanll2010HTasker fc Tan! 
2009: lTasker fc Brvanll2008L 120061) . 

We used a three-dimensional box of side 32kpc with 
a root grid of 256 3 and four levels of refinement, giv- 
ing a limiting resolution (smallest cell size) of 7.8 pc. 
For the effect of resolution on our simulation results, see 
TT09 where a detailed study is presented. Gas was re- 
fined whenever the Jeans' Length was r esolved by less 
than four cell widths, as suggested by iTruelove et al.l 
(|1997l) as the resolution needed for avoiding artificial 
fragmentation. On our finest level, the Truelove criteria 
is maintained until ~ 100 cm -3 , our threshold density for 
cloud definition. A discussion on resolving the gravita- 
tional collapse in the simulation is presented more fully 
inTT09. 

The evolution of the gas in Enzo was performed using 
a three-di mensional version of the Zeus hydrodynamics 
algorithm (St one fc Norm an 1993). This routine uses an 
artificial viscosity term to represent shocks, where the 
variable associated with this, the quadratic artificial vis- 
cosity, was set to 2.0 (the default) for all simulations. 

Radiative cooling w as included using rates from the 
analytical expression of lSarazin fc Wh ite (1987J) for solar 
meta llicity to 10 4 K and down to T = 300 K using rates 
from iRosen fc Bregmanl (| 19951 ) . This allows the gas to 
cool to the temperat ure of the upper en d of the atomic 
cold neutral medium (|Wolfire et alj|2003l ). Actual GMCs 
will have temperatures of ~ 10 K, an order of magnitude 
below our minimum temperature. However, at our reso- 
lution, clouds with diameters of 100 pc only have 13 cells 
in each linear dimension (with an average GMC having a 
radius of 16kpc and 4 cells across), which is insufficient 
to resolve the full turbulent structure of the gas. We also 
do not include pressure from magnetic fields, so impos- 
ing this temperature floor of 300 K produces a minimum 
sound speed of 1.8 kms -1 to crudely allow for these ef- 
fects. In fact, the velocity dispersion within our clouds is 
typically higher than this by about a factor of two, im- 
plying that this floor is not having a significant impact 
on our cloud properties. 

In addition to radiative cooling, the gas can also be 
heated via diffuse photoelectric heating in which elec- 
trons are ejected from dust grains via FUV photos. In 
the simulation where this was turned on, we included a 
radially dependent h eating term of the form described in 
iWolfire et all (l200l : 

_ ln -24 r / er {R ~ R °V HR ergss" 1 r > 4.0 kpc 

Ipe-lU £,,00 | e -(4-flo)/H Serg5S -l r<4 .0kpc 

(4) 



where the heating efficiency th — 0.05 and G n is the 
incide nt far-ultraviolet field normalized to the iHabind 
( 1968) estimate for the local IS M value . We t ake a value 
of Go = 1.7 in agreement with iDraind (|1978[ ). i?o is the 
radial scale length at 8. kpc and Hr = 4.1 kpc, the scale 
length as estimated by IWolfire et al.l (|2003l) . 

Collisionless star particles, representing star clusters, 
are allowed to form in our simulation in the main region 
of the disk between 2.5 < r < 8.5 kpc. As is described in 
12.21 below, this is the area of our disk where we identify 
and analyze the GMCs. Within this region, star parti- 
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cles are created when the density within a cell exceeds 
the threshold value of rva = 100cm -3 . Since our finest 
refinement cells (where the stars will be formed) are still 
7.8 pc across, the gas within them can be assumed to be 
turbulent. We therefore do not check for gravitational 
collapse or boundedness of the cell gas, since such pro- 
cesses are likely to affect star formation on much smaller 
scales than what we can resolve. Cells with tempera- 
tures greater than 3000 K are also prevented from form- 
ing stars to rule out the possibility of star formation in 
the hot dense gas of shock fronts. In practice, star par- 
ticles form in gas typically close to 300 K, so this limit 
is not particularly important. When a cell reaches the 
threshold density, a star particle is created whose mass 
is calculated by: 



Iff 



(5) 



where e is the star formation efficiency (the fraction of 
gas that is converted into star particles per dynamical 
time), At is the size of the time step, ts is the time for 

dynamical collapse in the cell (tg = ( J ) and p gas 

is the gas density. The resultant object should be consid- 
ered a star cluster, rather than an individual star since 
it contains approximately 1000 M Q . For our simulations, 
we chose a value for eg = 0.02, in ag reement with the ob- 
servat ional constraints described bvlZuckerman fc Evaiisl 
(1971 for GMCs and iKrumholz fc McKed (|2005t ) for 
GMCs and their internal, higher density components. 

An additional computational requirement is that the 
star particle will not be created if the calculated value for 
is less than a given minimum value of m m i„ = 10 3 M Q . 
This purely numerical addition is to prevent the calcu- 
lation becoming prohibitively slow due to a extremely 
large number of star particles. In the situation where 
this occurs, an override exists whereby a particle of mass 
w m i n is created with a probability equal to the ratio be- 
tween the mass of the would-be star particle and m m i n - 
In practice, almost all star particles born in the high- 
est resolution simulations are created via this stochastic 
method with masses equal to m m in- 

The motions of the star particles are calculated as a col- 
lisionless N-body system. They interact gravitationally 
with the gas via a cloud-in-cell mapping of their positions 
onto the grid to produce a discretized density field. The 
number of star particles created during the simulations 
is 2.5 - 3 million. In this paper, there is no localized 
feedback to the gas directly from the star particle. 

2.1.1. Runs Performed 

This paper presents the results from two different sim- 
ulations, both at a limiting resolution of 7.8 pc. 

The first simulation (simulation disk SFOnly) includes 
star formation, implemented as described above, but 
with no form of additional heating. The second sim- 
ulation (simulation disk SF+PEheat) also includes star 
formation and a radially dependent diffuse heating term 
with the form given by Equation |4] 

The results in the absence of both star formation and 
diffuse heating are presented in TT09, simulation disk 
NoSF. 



2.2. Galaxy Initial Conditions and Cloud Analysis 

The initial conditions for the simulations are described 
in detail in TT09. They consist of an isolated gas disk 
sitting in a static background potential that represents 
both a dark matter halo and a stellar disk component for 
a galaxy similar to the Milky Way. The potential gives 
the disk a constant circular velocity for r » 0.5 kpc of 
v c = 200 km s" 1 . 

We focus on the gas in the main region of the disk be- 
tween 2.0 < r < 10.0 kpc. Gas here is initially marginally 
stable against gravitational collapse, wit h the Toomre 
Q parameter for gravitational instability (|Toomrelll964l) 
having a constant value: 



ttGXL 
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(6) 



where k is the epicycle frequency and E a , the gas sur- 
face density. a g = y/ of lt + c 2 s is the mass- weighted ID 
velocity dispersion of the gas, with cr nt being the veloc- 
ity motions in the disk plane after subtraction of the 
circular velocity. The vertical profile is proportional to 
sech (z/z/s), where the scale height, Zh varies with galac- 
tocen tric radius based on HI ob servations of the Milky 
Way (Bin nev k. Merrifierdi n~998). At the solar radius of 
8 kpc, Zh = 290 pc. For a flat rotation curve, k = \/2v c /r, 
which gives a gas density profile of the form: 



p(r,z) 
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Low density regions of gravitationally stable gas sit in 
the disk center between < r < 2 kpc and at the outer 
edge between 10 < r < 12 kpc. 

As the gas cools, the main region of the disk becomes 
gravitationally unstable with Q < 1 and fragments. Star 
formation and cloud analysis are restricted to between 
2.5 > r > 8.5 kpc; within the main region but avoiding 
boundary effects. 

Details of the algorithm used to identify and track 
the GMCs in the simulations are described in TT09. 
In short, we identify "GMCs", i.e. star-forming clouds, 
as peaked and coherent structures contained within con- 
tours of the threshold density of nn,c > 100 cm -3 , about 
the mean volume density of observed galactic GMCs. It 
is worth remembering that while GMCs are by definition 
molecular, our gas is purely atomic, so we are selecting 
structures that would be mostly molecular in reality, al- 
though our procedure does not distinguish dense atomic 
gas that might be present in photodissociation regions. 
By comparing outputs of the simulation at 1 Myr inter- 
vals, the clouds are tracked over the course of the simu- 
lation to produce a timeline of their evolution. The sim- 
ulations themselves are run for 300 Myr, just over one 
orbital period at the outer edge of the main region in the 
disk. 

3. GLOBAL EVOLUTION OF THE GALACTIC DISK 

The disk rapidly cools from its initial conditions to 
fragment through gravitational instabilities (see TT09 
for a discussion of this process). By t = 140 Myr, the 
main region of the disk between 2.5 > r > 8.5 kpc has 



Fig. 1. — Images of the galactic disk at 200 Myr. Top panel is for disk SFOnly while the bottom panel shows disk SF+PEheat. The 
left-hand image is 20kpc across and shows the log-scaled surface density of the disk with range [0.0033,3715] Mqpc~ 2 . The center pane 
is the projected star particle density and the right-hand image shows a 2 kpc density slice of the mid-plane. Blue contour lines mark the 
cloud boundaries corresponding to a number density of 100 cm -3 and the new star particles (age < 1 Myr) are shown in red. 



fully fragmented into dense cold clouds of gas, embedded 
in a warmer medium. 

Three images of the disk at t = 200 Myr are shown in 
Figure Q] for simulation disk SFOnly (top row) and disk 
SF+PEheat (bottom row), where diffuse heating is in- 
cluded. The left-hand panel shows the gas surface density 
over the main region of the disk, the middle panel shows 
the star particle density for the same region and the 
right-hand panel displays a 2 kpc slice of a typical patch 
in the galactic mid-plane. In this right-hand panel, blue 
contours mark our cloud threshold definition of 100 cm" 3 
and new star particles, with age < lMyr, are shown in 
red. 

The effect of including the diffuse heating on the disk's 
ISM can be seen in the images of the gas surface density. 
Without the additional heating in disk SFOnly, the gas 
collapses to form smaller structures, especially in the in- 
ner region, 2 > r > 4, where both the density of the disk 
is highest (and hence the dynamical time shortest) and 
the heating term, as described in equation 0] is at its 
maximum. While the disk has evidently become grav- 
itationally unstable in both cases, its fragmentation is 
reduced by the diffuse heating. 

One of the main results of the reduced fragmentation 
is shown in the middle panel of Figure [T] in the distribu- 
tion of star particles. The stellar density is visibly lower 
when diffuse heating is included. At t = 200 Myr, disk 



SFOnly has formed over 3 million star particles, whereas 
this is reduced to 2.6 million in disk SF+PEheat. The 
star formation properties of the disk will be considered 
quantitatively in section [6] 

The findings from the first two image panes in Fig- 
ure Q] are supported in the close-up of the 2 kpc slice of 
the galactic plane on the far-right. The filaments con- 
necting sites of star formation are denser and thicker in 
disk SF+PEheat and the smaller star formation regions 
appear to have less stars within them than ones of similar 
diameter in disk SFOnly. 

The new star particles are born inside the identified 
cloud boundaries, as expected since the density thresh- 
old for both cloud identification and star formation is 
njj.c = 100 cm~ 3 . In the occasional case where a star 
particle appears to exist just outside the cloud boundary, 
it has either moved during the last 1 Myr or destroyed 
(via mass removal) the small part of its cloud that it was 
formed in. With diffuse heating, the gas between the 
clouds forms thicker filaments, supported against further 
fragmentation by the thermal addition. This gives the 
warm ISM between the clouds a more coherent struc- 
ture. 

This variation in the global structure of the disk, to- 
gether with its evolution, is shown quantitatively in Fig- 
ure [2] The four plots show azimuthally averaged radial 
profiles of the disk properties for disk SFOnly (top row) 
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Fig. 2. — Evolution of the azimuthally-averaged radial galactic profiles for disk SFOnly (top row) and disk SF+PEheat which (bottom 
row). Plots left to right show: (a) gas mass surface density (S 9 ), (b) ID velocity dispersion of the gas (c 9 ), (c) temperature (T) and (d) 

Toomre Q parameter. Note, S 9 = f^^pc P( z )d- Z i T is a mass-weighted average over —1 kpc < z < 1 kpc and Q makes use of <j g evaluated 
as a mass-weighted average over the same volume. 



and disk SF+PEheat (bottom). From left to right, the 

plots are: (1) gas surface density, S g = f^ lk ^ p(z)dz, 

(2) the gas velocity dispersion, <r s , calculated via a mass- 
weighted average over — 1 kpc < z < 1 kpc utilizing only 
disk plane velocity components, (3) the mass-weighted 
temperature, T, averaged over the same range as a g and 
(4) the Toomre Q parameter for gravitational instabil- 
ity as given in Equation |6] The four lines on each plot 
show the profile at different times (t = 0, 100, 200 and 
300 Myr) during the simulation. 

The initial conditions at t = are shown by the black 
solid line in all plots. We can see our main disk region 
between 2.5 < r < 8.5 kpc is initially at a surface density 
higher than the surrounding gas by over an order of mag- 
nitude and has a constant borderline stable Q — 1.5. The 
initial temperature is constant over the entire simulation 
box and corresponds to a sound speed of c s = 9kms _1 . 

During the first 100 Myr (red dashed line) the gas cools, 
bringing Q below the critical stability value of 1.0 and 
causing the gas to fragment. In the top panel showing 
disk SFOnly, the temperature drops until it approaches 
the floor of our radiative cooling curve at 300 K, whereas 
the addition of diffuse heating in the bottom panel slows 
down the cooling of the densest clumps, increases the av- 
erage minimum temperature reached by approximately a 
factor of 2.0. Global ring instabilities appear in the sur- 
face density profiles of both simulations (see TT09 for 
images of their formation) , but are less prominent in the 
gas warmed by the diffuse heating, which provides an ex- 
tra thermal support. As the disk fragments tangentially, 
these fluctuations flatten out. The magnitude of the sur- 
face density remains almost constant over the disk's main 
region during this early period. At lower radii, how- 



ever, the circular motion becomes poorly resolved by the 
Cartesian grid, causing an infall of low density gas from 
the inner region of the disk to build up at the disk center. 

By the time the simulation reaches 200 Myr (blue dot- 
dashed line), star formation has depleted the gas, caus- 
ing the surface density to drop by a factor of 2.0 in disk 
SFOnly and 1.5 for disk SF+PEheat. Since star for- 
mation occurs in the coldest and densest regions of the 
disk, the azimuthally averaged temperature rises as do 
the corresponding values for the velocity dispersion and 
Q. The velocity dispersion also increases due to gravita- 
tional heating, as star particles are formed to create local 
deep potential wells. This effect is strongest in the disk 
without diffuse heating, since 500,000 more star particles 
have formed by this time, causing a greater depletion of 
the dense gas. 

Star formation continues to deplete the gas, resulting 
in the surface density dropping by a factor of about 6 
from the initial conditions by 300 Myr at r = 6 kpc for 
both runs (green dotted line) . The SFR is highest in the 
densest gas (since it is proportional to the free-fall time) , 
causing the surface density to drop first in the inner part 
of our main region, 2.5 < r < 4 kpc. This is marked most 
in disk SF+PEheat, where the star formation progresses 
more slowly due to the reduced fragmentation. 

The velocity dispersion is significantly lower in disk 
SF+PEheat, having increased more slowly over the 
300 Myr. This initially seems surprising, since the tem- 
perature is higher, but diffuse heating has reduced the 
ability of the disk to fragment into massive bound clus- 
ters, as was seen in Figure [TJ leaving a filamentary state 
that reduces the in-plane velocity motions, cr n t, (and the 
corresponding Q value) , that are excited by gravitational 
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interactions between fragmented clumps. 

4. THE STRUCTURE OF THE ISM 

A one dimensional analysis of the structure of the ISM 
in the main region of the disk can be seen in the proba- 
bility distribution functions (PDFs) plotted in Figure |3] 
for the simulations at the same times shown in Figure [2l 
The top two graphs show the gas volume fraction as a 
function of density for the disk SFOnly (left) and for disk 
SF+PEheat (right). The bottom graphs show the mass 
fraction over the same density range. 

Without diffuse heating, the mass fraction of gas in 
disk SFOnly above our cloud definition limit of nu. c > 
100 cm -3 is 0.43, 0.4 and 0.16 for simulation times t = 
100, 200 and 300 Myr respectively. For disk SF+PEheat, 
the same output times have cloud mass fractions of 0.36, 
0.46, 0.33. The rise in dense gas from 100 to 200 Myr 
when diffuse heating is present is evidence for the disk 
still fragmenting over this period. 

Little evolution in the shape of the PDFs is seen over 
the course of the simulations once the initial condi- 
tions have produced a fragmented disk. In contrast to 
disk NoSF in TT09, we do not see a rise in high den- 
sity gas over time, since stars now form in gas above 
Uh,c > 100 cm -3 and deplete the gas reservoir above this 
threshold. This will have an impact on the maximum 
cloud mass as will be seen in section [5j 

In agreement with other simulations of galactic disks 
(e.g.lRobertson fe Kravtsovll2008t iTasker fe Brvanll2008t 
iWada fe Normanll2007D . both our disks can be fitted with 
a log-normal tail to their volume- weighted PDFs. In disk 
NoSF in TT09, the fit to the 200 Myr and 300 Myr lines 
extended up to densities of rin < 10 5 cm -3 , but in these 
cases star formation removes gas to steepen the profile 
at the high density tip, making this fit not as good. This 
log normal fit from disk NoSF is shown in Figure [3] and 
has the form, 

PDF = _ e -(lna;-lna;) 2 /2gg DF ^ 

(7pDFV27r 

where x — p/p and upbf = 2.0. It remains a good fit 
to both disks at t = 200 Myr, but by t = 300 Myr, star 
formation has significantly eroded the dense gas. 

The one-dimensional Mach number, Ai, for the star- 
forming high density tail of the PDF can be estimated 
from the azimuthally averages profiles in Figure [2J Be- 
fore star formation significantly depletes the disk gas, the 
velocity dispersion, a g , is dominated by the in-place ve- 
locity motions, cr n t, which are of order 15kms _1 . In this 
high density region, the temperature of the gas is close to 
our cooling floor at 300 K (a fact that will also be seen in 
Figured]), giving a c s = 1.8 kms -1 . The Mach number is 
given by the ratio of these two values, M = er nt /c s « 8.3. 

From this and the value of (Jpdf, the nature of the 
turbulence production in the star-forming gas can be 
deduced. The two numbers are related via Cp DF = 
In [l + b 2 A4 2 ] , where b is found to vary between b ~ 
1 /3 for solenoidal (divergence- free) turb ulent modes and 
6 ~ 1 for compressive (curl-free) modes (F ederrath et al.l 
120081) . For M = 8.3 and a PDF = 2.0, we find b = 0.88, 
suggesting that compressive forcing is the dominant tur- 
bulent mode. This is consistent with the compressive 



nature of cloud collisions which are likely to be a driving 
force for the turbulence. 

Removing gas from the high density tail of the PDF 
via star formation causes an increase in the fraction of 
gas at densities below nn < lcm -3 . This effect is more 
pronounced in disk SFOnly, since the fraction of cloud 
gas has dropped by almost a factor of 3 between 100 and 
300 Myr, compared with a maximum change of just 1.4 
in disk SF+PEheat. 

Despite its temporal continuity, the fraction of gas be- 
tween 10~ 4 < rin < 1 cm~ 3 in disk SF+PEheat is higher 
than for disk SFOnly, remaining close to the initial condi- 
tion value. This gas is the warm ISM in which the clouds 
are embedded and, as seen in Figure [TJ has a denser fil- 
amentary structure due to a smaller level of fragmenta- 
tion. 

The peak at low densities of nn < 10~ 4 cm~ 3 corre- 
sponds to gas in the 1 kpc region above and below the 
disk plane. In the disk SFOnly, we see the density of 
this region increase over time as the high velocity dis- 
persion seen in Figure [5] causes the disk scale height to 
increase. With diffuse heating, the disk settles to an ini- 
tially thicker profile with a lower velocity dispersion in 
its less fragmented ISM, undergoing less evolution. 

A two-dimensional representation of the evolution of 
the ISM is shown as density vs temperature contour 
plots in Figure |U The upper two rows of plots show 
the gas volume distribution (top) and gas mass distribu- 
tion (bottom) at times t = 100, 200 and 300 Myr for disk 
SFOnly while the lower six plots show the same distri- 
butions for disk SF+PEheat. Diagonal lines mark con- 
stant pressu re with the values bein g representative of the 
Milky Way (|Boulares fe Coxlll990t ): The upper-most line 
is the estimated value for the total Milky Way pressure, 
Ptot/kb = 2.8 x 10 4 K cm~ 3 . The lower two lines show 
the thermal pressure, Pth/kb — 0.36 x 10 4 K cm" 3 , and 
the thermal pressure excluding the hot gas component, 
Pth,„ohotAh = 0.14 x 10 4 K cm" 3 . 

The distributions show that gas warmer than 10 4 K, 
the typical temperature fo r the warm interstellar medium 
(jMcKee fc Ostrikerlll977[ ). lies largely in pressure equi- 
librium, although with a wide range of values that run 
over a continuous distribution of densities and tempera- 
tures. Without diffuse heating, gas in disk SFOnly that 
is cooler than this remains in pressure equilibrium until 
it approaches the floor in our radiative cooling curve at 
300 K. This appears as a sharp spike in high mass and low 
volume density in all plots. The gas in this region is in 
the star forming clouds whose self-gravity causes them to 
be over pressurized with respect to the surrounding ISM. 
There is a perceptible decrease of gas in this spike over 
time as the cloud gas is converted into stars. Tempera- 
tures lower than the cooling curve minimum are achieved 
via adiabatic expansion. 

With the addition of diffuse heating in disk 
SF+PEheat, cooler gas is warmed causing its pressure 
to increase. This results in a significantly thinner profile 
and gas below T < 10 4 K being out of pressure equilib- 
rium. The smaller range of temperatures and densities 
in the warm ISM is indicative of the coherent structure 
of filaments we saw in Figure [1] and the lower velocity 
dispersion in Figure [2j 

Overtime, both disk profiles broaden as cloud mergers 
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Fig. 3. — Probability distribution function (PDF) for the disks over the radii 2.5 < r < 8.5 kpc and height —1 kpc < z < 1 kpc. 
Left-hand plots are for the disk with only star formation (disk SFOnly) right-hand plots are for the disk where diffuse heating is included 
(disk SF+PEheat). The top panel shows the evolution of the volume- weighted PDF where the solid- line curve shows a log normal fit to 
the high density PDF tail. The bottom panel is the mass-weighted PDF. There is relatively little evolution in the PDF shape over the 
course of the simulation, although the depletion of gas in the SFOnly disk is evident in the gas fraction shift to lower densities. This is not 
as evident in SF+PEheat and more mass exists in the mid-density (l0 -4 < nu < 1cm -3 ) ISM. 



and tidal encounters disrupt the structure of the dense 
gas. 

For the gas in pressure equilibrium at T > 10 4 K, 
its value is lower than that of the Milky Way in both 
the disks during the majority of the simulation. This 
is not surprising, since gas at this temperature will be 
strongly affected by sources of energetic feedback which 
we do not consider. We do see a rise in the pressure 
over time, resulting in a factor of 10 increase between 
t = 100 to 300 Myr. Comparing with the results in disk 
NoSF, where no such increase was seen, we conclude this 
is the result of the localized potential created by the star 
particles. 



5. THE PROPERTIES OF THE CLOUDS 

While the galaxy's global structure concerns the en- 
tire ISM, its star formation occurs almost solely in the 
densest, coldest component of that gas. This cold phase 
is organized into the extended structures of the GMCs. 
As the disk cools and fragments, regions of gas exceed 
our threshold of 100 cm -3 and we recognize them as the 
GMCs. The properties of these entities dictate the envi- 
ronment in which stars will form and their birth, evolu- 
tion and death will determine the star formation rate in 
the galaxy. 

In this section, we focus on the properties of the iden- 
tified GMCs, examining the the evolution of the clouds 
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4. — Density vs temperature contour plots for disk SFOnly (upper panel of six images) and disk SF+PEheat (lower panel of six), 
panel, the top row shows the distribution of gas volume while the bottom row shows gas mass for the region 2.5 kpc < r < 8.5 kpc 
kpc < z < 1 kpc. Solid lines show the total pressure in the Milky Way, Ptot/^b = 2.8 X 10 4 K cm -3 (top), the thermal pressure, 
= 0.36 X 10 4 K cm 3 (middle) and the thermal pressure excluding the hot gas component, P t ^ no hot/fcf) = 0.14 X 10 4 K cm 3 
i) lIBoulares fc Coxlll"990T) . We see a continuous range of densities and pressures in both disks that is in pressure equilibrium above 
itures of 10 4 K and over-pressurized in the self-gravitating clouds. Diffuse heating greatly reduces the range of values due to the 
^locity dispersion and denser warm ISM. 
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Fig. 5. — Left: Evolution of the total number of GMCs, N c , (on left axis) and total number of star particles, N stlirs , (right axis) for both 
disks. Right: the GMC formation rate, dN c j orm /dt, and the GMC merger rate (equivalently, the rate of cloud destruction via mergers, 
since these are virtually all between two clouds), |<i/V Cimcr g C / 'dt\ for both disks. The disk fragments over the first 100 Myr where the number 
of clouds begins to decrease due to mergers between clouds and gas depletion from star formation. The latter process is slower in the gas 
with diffuse heating, which suppresses the fragmentation of the disk and star formation in the clouds.) 
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Fig. 6. — Cloud merger timescales, averaging over 50 Myr intervals of simulation time, for disks SFOnly (left) and SF+PEheat (right). 
Only clouds born after 140 Myr, the initial fragmentation of the disk, are included in the analysis. The average time for a merger is low 
at ~ 0.25 of an orbital period and largely independent of galactocentric radius. The dashed line shows the merger times of clouds with 
M c > 10 6 Mq (i.e., the average time for a cloud of that size to undergo a merger with a cloud of any mass), evaluated over the interval t 
= 225-275 Myr, which is lower by a factor of 2. 



11 



both as a function of simulation time, i, and as their age, 
or cloud time, t' . 

The evolution of the cloud population itself is shown 
Figure [5] for both disks. The left-hand plot shows the 
total number of clouds and stars in the disk, while the 
right-hand plot shows the formation and merger rate of 
the clouds. The number of clouds initially increases as 
the disk gravitationally fragments, peaking at a rate of 
around 400 Myr -1 at 100 Myr in both disks. After that 
time, the disk is fully fragmented and both the cloud 
population and formation rate drop, accompanied by an 
increase in the number of star particles, as gas is con- 
sumed during star formation and mergers further erode 
the cloud population. There is a delay in the peak of 
the cloud formation rate and cloud merger rate that is 
equivalent to the mean merger rate between clouds. By 
300 Myr, the cloud population has dropped by a factor of 
10 without diffuse heating and a factor of 8 where diffuse 
heating was included. When diffuse heating is present, 
the build up of the stellar population is slightly more 
gradual, indicative of a slower star formation rate. By 
300 Myr, the total number of star particles in both disks 
is approximately the same at ~ 4 million, but at 200 Myr, 
it differs by ~ 20%. This accounts for the higher number 
of clouds in disk SF+PEheat at 300 Myr; the slower star 
formation rate allows the cloud time to accrete material 
as it is consumed, allowing a greater number to survive 
despite the ultimate star particle number being the same 
as in disk SFOnly. 

Mergers between clouds happen throughout the course 
of the simulation, their rate related to the number of 
clouds in the disk. How mergers are recorded by the 
tracking algorithm is discussed in TT09, but in short, a 
merger is said to have happened when a single cloud is 
at the predicted position for two other clouds after 1 Myr 
of evolution, within a margin of twice the average radius 
of the potential merger product. 

In the disk SFOnly, the number of mergers drops more 
steeply in the second half of the simulation compared 
to SF+PEheat. This corresponds to a steeper fall off 
in the number of clouds as the mass of gas above the 
cloud threshold density drops heavily as was seen in Fig- 
ure [3l As we will see in Section |6l gas depletion in this 
simulation also causes the star formation rate to drop 
considerably over this period compared to the diffusely 
heated disk. 

The frequency of the mergers in both disks are shown 
in Figure [6] as a function of the galactocentric radius. 
The merger time is computed as a fraction of the orbital 
time and averaged over 50 Myr. Merger rates over four 
time intervals during the simulation (75-125 Myr, 175- 
225 Myr, 225-275 Myr, 275-325 Myr) are plotted, where 
a merger is considered to be between two clouds. Cloud 
mergers involving more than two objects are negligibly 
rare. 

There is no trend between merger rate and galacto- 
centric radius for the majority of the simulation. At late 
times, disk SF+PEheat undergoes less mergers in its in- 
ner 4kpc, corresponding to the drop in gas density (and 
so cloud number density) in that region, coupled with 
the lower velocity dispersion, as seen in Figure [5] In disk 
SFOnly, there is a spatially uniform increase at latest 
simulation time in the time between mergers, due to the 
reduced number of clouds seen in Figure [5j 



Within the first three time frames, the average num- 
ber of mergers is approximately constant in both simula- 
tions, with a value of ~ 0.25 of the orbital time, close to 
that seen in TT09 for disk NoSF. For our largest clouds, 
M > 10 6 M Q , this rate is lower by a factor of 2 . This is in 
agreement with the calculations of lTanl (|2000h . who esti- 
mated frmergerAorbit = 0-2 in an analytical model. Since 
the disk circular velocity is v c = 200 kms" 1 , the average 
time between interactions is ~ 30 Myr at a radii of 4kpc. 
This is slightly higher than for the clouds in disk NoSF 
where an average time of 25 Myr was found, due to the 
star formation depleting the number of clouds. 

The frequency of the cloud collisions is indicative that 
they play a major role in shaping the properties of the 
GMCs. What is less clear is what effect they have on the 
star formation rate as will be discussed in Section [6] 

5.1. GMC Properties with Simulation Time 
5.1.1. Simulation Time: Cloud Property Distributions 

The properties of the individual clouds are shown in 
Figure [7] for disk SFOnly and Figure [8] for the clouds for 
disk SF+PEheat. The distributions are shown at three 
different times during the simulation; 100 (red dashed 
lines), 200 (blue dot-dashed lines) and 300 Myr (green 
dotted lines). 

The top left plot, Figure [3(a) and Figure [8](a) , show 
the profiles for the cloud gas mass. In difference to when 
no star formation occurred in disk NoSF and in agree- 
ment with the PDF plots in Figure El we do not see an 
increase in the cloud mass over the course of the sim- 
ulation. Instead, star formation prevents the formation 
of a high mass tail and caps the maximum cloud mass 
at M < 1O 6,7 M . When diffuse heating is included 
in Figure [5](a), we see one outlier beyond this trunca- 
tion at M = 10 71 M Q which could be due to a recent 
merger event. Otherwise, the addition of heating does 
not change the mass cut-off. 

The peak of the mass distribution appears at 10 5 ' 6 M 
at 100 Myr after the start of the simulation in disk 
SFOnly and slightly lower, 10 5 5 M , in disk SF+PEheat. 
This drop in peak mass is due to the heating term sup- 
porting the cloud against its own gravity, slowly the col- 
lapse to its maximum density. This was noted in the 
disk PDFs in Figure [3] and will be seen again below in 
the distribution of cloud surface densities. 

Without diffuse heating, the conversion of gas into 
stars causes the peak to migrate to lower masses over 
time, moving to 10 49 M Q by 300 Myr. At this stage, 
75 % of the clouds have a stellar mass fraction greater 
than 90%. This high stellar content is due to the lack 
of feedback that would act to prevent gas continually 
collapsing to form stars within the cloud. With diffuse 
heating, the evolution of the mass profile is reduced, with 
no decrease in the peak position over time and a smaller 
evolution in the low mass end of the spectrum. This is in 
keeping with Figure [3] where we saw little change in the 
heated disk's ISM over time. Again, the added thermal 
support is the cause here, decreasing the cloud density 
and so increasing its dynamical time, causing a slower 
rate of conversion from gas into stars through equation[5l 

Obs ervations of the Milky Way by Wi lliams fc McKeel 
( 1997) find a truncation in the cloud molecular mass 
at M u ~ 6 x 1O 6 M . Since we follow all the neutral 
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Fig. 7. — Normalized distributions of GMC properties at 100, 200 and 300Myrs after the start of the simulation for disk SFOnly. 
Top left plot shows: (a) cloud mass distribution overlaid with fits to observational data from the Milky Way (blue and violet solid lines; 
IWilliam s & McKcc (1997)). Middle left, (b): the average radii of clouds, R c a = (^c/t) 1 / 2 where A c is the projected area of the clouds 
in the Y-Z plane. Bottom left, (c): mass surface density of the clouds, S c j\. Top right, (d): virial parameter, a vlr with the vertical line 
indicating where ce vlr = 1, the limit for gravitational binding. Middle right, (e): vertical component of the specific angular momentum 
vector, j' z . The vertical line indicates a value of j z of a spherical (~ 110 pc radius) region of the initial conditions at galactocentric radius 
r = 4kpc containing 10 6 Mq. Bottom right, (f): vertical height distribution, z, of the clouds. 
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Fig. 8. — Normalized distributions of GMC properties at 100, 200 and 300Myrs after the start of the simulation for disk SF+PEheat. 
Plots are the same as those described in Figure[7] Top left: (a) cloud mass, M c . Middle left: (b) cloud radius, Rc,A- Bottom left: (c) 
mass surface density, S c . Top right (d) virial parameter, o v ir- Middle right (e) vertical component of the specific angular momentum, j z . 
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gas, our results also i nclude an atomic envelope whose 
mass is estimated by iBlitz et all (|1990l ) to be of order 
the molecular component. This would take the observed 
truncation li mit up to M < 1.2 x 10 7 M Q . Recent ob- 
servations by iFukui et all (|2009h of the LMC, suggest 
smaller envelopes with 10% of the mass, providing a limit 
of M u ~ 6.6 x 10 6 M Q . In both cases, our maximum mass 
falls below or close to the observed truncation. 

This simulated mass spectrum can be compared to ob- 
servations of GMCs in the Milky Way. These observa- 
tions (and i ndeed, those in other ga laxies such as M33 
as found by IRosolowskv et afl [2003) fit a power law of 
the form oc M c Qc . The two lines on Figure 0(a) 

and Figure [5](a) show this re lation for observatio n s of a c 
in the Milky Way made by IWilliams k McKed (|1997t ). 
Their measured value for a c depends on if the observed 
sample is assumed to be equally under-sampled at all 
masses (a c = 0.6, MWA) or if lower mass clouds are 
more under-sampled than high mass (a c = 0.8, MWB). 

The distributions most appropriate to compare with 
these observations are the lines at 200 and 300 Myr, after 
the disk has settled and cloud interactions have become 
the dominant force in the ISM. Within the mass range 
~ (0.5 — 10) x 10 6 Mq, where our clouds are most re- 
solved, we find that disk SFOnly approaches the Milky 
Way GMC population. With diffuse heating, the evolu- 
tion of the mass profile is reduced, and the distribution 
is steeper at later tim es. 

IWilliams fc McKed (|1997l ) estimate that, within the 
inner Milky Way, there are ~ 1000 clouds with M c > 
10 5 M o and 100-200 clouds with M c > 10 6 M o . At 
200 Myr, we find for disk SFOnly 1122 clouds with 
M c > 10 5 M Q and 179 clouds with M c > 10 6 M o , in 
good agreement with this result. The cloud population 
in disk SF+PEheat is slightly larger with 1491 clouds 
with M c > 10 5 M Q and 238 with M c > 10 6 M Q , reflect- 
ing the higher fraction of cloud gas at this time, but is 
still of the same order as the observations. Moreover, 
the mass fraction of clouds at 200 Myr is 0.4 for the disk 
without diffuse heating and 0.46 when it is included (as 
se en in section 4 |, coi nciding with the estimated value 
of IWolfire et al.l (|2003|) for the fraction of molecular to 
atomic gas inside the solar radius. 

Figure 0(b) and Figure [8](b) show the average radius 
of the clouds defined as R c ,a = {Ac/it) 1 / 2 , where A c is 
the projected area of the cloud in the y — z plane (that is, 
as it would be measured by an observer embedded in the 
plane of the galaxy). Cloud radii are typically around 
15 pc in both disks, in good agreement with observations 
of GMCs in the Milky Way and M33 (|Lada et al.ll201fl 
IRosolowskv et al.l 12003) . This is slightly smaller than 
the non-star forming clouds in disk NoSF (typical radius 
~ 20 pc), but the inclusion of star formation has only sig- 
nificantly impacted the tail of the distribution which, as 
with the mass, is now prevented from continual growth. 

As with the mass profile in Figure 0(a), there is a de- 
crease in the peak cloud radius over time for disk SFOnly, 
as gas is converted into stars. This is not reflected in disk 
SF+PEheat as its cloud mass fraction changes slowly un- 
der the influence of a lower SFR. At later times, there are 
a small number of extended structures which are likely 
to be from older clouds that have undergone multiple 
merger events. 



Figure0(c) and Figure[8](c) show the surface density of 
the clouds, S c = M c /A c . The distribution is the most ro- 
bust property plotted, showing minimum evolution over 
time. The peak value sits ~ 300M©pc -2 , in agreement 
with disk NoSF, although at 100 Myr, disk SFOnly has a 
lower higher peak surface density than disk SF+PEheat, 
as the heating slows down the fragmentation to the high- 
est densities. 

There is a minor amount of evolution over time. With- 
out diffuse heating, the surface densities drop as the 
clouds becomes more stellar dominated. With the diffuse 
heating, the shift is in the opposite direction as the disk 
self-gravity slowly overwhelms the support from heat- 
ing and amount of cloud gas rises between 100 Myr and 
200 Myr. Compared to the changes in mass and radius 
profile, however, these shifts are small. 

The independence of the surf ace density p rofile 
to physics an d time agrees with iLarsonl (|1981l) and 
ISolomon etafl ((1987f t. who found that GMCs in the 
Milky Way had a constant surface density. Their mea- 
surements give a value for E c lower than what we find, 
with S c = 200M©pc~ 2 , but this does not take into ac- 
count the atomic envelope, which c an be a factor of 100 
lower in density (|Fukui et al.ll2009D . 

The degree of gravitational binding in the clouds can be 
estimated by the alpha virial parameter, a rmv i r , plotted 
in Figure0(d) and Figure[8](d). This quantity is defined 
as a V i r = 5(T 2 R Ct A/(GM t ), where a c is the mass averaged 
velocity dispersion of the cloud, a c = (c^ + cr 2 ^) 1 / 2 , with 
Cnt,c the one-dimensional rms velocit y dispersion about 
the c loud's center-of-mass velocity ([Bertoldi fc McKeel 
I1992D . In contrast to TT09, M t is now the sum of the 
total gas and stellar mass in the cloud. To calculate this, 
star particles are associated with the cloud if they reside 
within its boundaries or within a distance of twice its 
average radius. A value of a V ir = 1.0 states that a spher- 
ical, uniform cloud with negligible surface pressure and 
magnetic fields is virialized. This translates to a cloud 
for which a v ir < 1-0 being gravitationally bound and one 
for which a v i r > 1-0 likely to be unbound. Observation- 
ally, clouds are seem to be on the borderline o f the two 
states, with a vir ~ 1 (jMcKee fc Ostrikerll2007l ) . 

For both our disks, the profile peaks at log(a^i r ) ~ 
—0.25 throughout the simulation, implying our clouds 
are weakly bound. There is an increase in the fraction 
of more highly bound objects over time as the clouds 
become more stellar dominated, resulting in mergers in- 
creasing the cloud's total mass more than its gas fraction 
(which would increase its radius or surface density) . This 
effect is strongest in disk SFOnly, with the distribution at 
300 Myr showing an almost flat profile below a v ; r < 0.5. 
As previously mentioned, 75% of the clouds at this time 
are now over 90% stellar in mass, the relatively small 
quantity of gas making it unlikely that our expression 
for a v i r can still apply. With diffuse heating this frac- 
tion is only ~ 46% and 300 Myr, and a v ; r shows a smaller 
level of evolution. 

The specific angular momentum in the vertical, z, com- 
ponent is shown in Figure 0(e) and Figure [11(e). The 
vertical line marks the specific angular momentum that 
a sphere of 10 6 M Q and radius r = 4kpc has in the ini- 
tial conditions. The clouds have an angular momentum 
much smaller than this, since they are much more com- 
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pact objects and therefore less susceptible to the disk 
sheer. 

Without diffuse heating, the angular momentum pro- 
file in disk SFOnly is very similar to disk NoSF, with 
j z decreasing over time, moving from an almost purely 
positive valued distribution to a broader profile where 
a substantial fraction of the clouds have a negative j z . 
These clouds rotate in the opposite sense to the galaxy as 
will be discussed in section ["5.1.3l This broadening is due 
to interactions between clouds affecting their rotation. 

The diffuse heating once again moderates the evolution 
of the cloud distributions, causing j z in disk SF+PEheat 
to show a smaller change over time. While cloud interac- 
tions are still common (as we saw in Figurc[5|) , their effect 
on the cloud's rotation is moderated by the denser, fila- 
mentary structure of the surrounding ISM which main- 
tains a low velocity dispersion as was seen in Figur^l 
This allows late forming clouds to form in a more struc- 
tured, less turbulent environment than without diffuse 
heating, which plays a significant role in determining the 
rotation of the clouds. 

The vertical distribution of the clouds is shown in the 
final panel, Figure [3(f) and Figure [8](f). The cloud scale 
height grows over the course of the simulation in both 
disks, due to the frequency of cloud collisions scattering 
clouds out of the galactic mid-plane. With the excep- 
tion of our last time analysis which has a high degree 
of scatter, the scale height of the clouds is compara- 
ble to that of the GMCs in the Milky Way at < 35 pc 
(|Stark &: Led 12005). Diffuse heating makes little differ- 
ence to this distribution, although we saw in Figure 
that the scale-height of the warm ISM is initially larger 
in disk SF+PEheat due to the added pressure from the 
heating. 

5.1.2. Simulation Time: Velocity Dispersion vs Size 
Relation 

The internal cloud velocity dispersion plotted against 
cloud size is shown for clouds present at t=200Myr in 
Figure |9] for both disks. Only clouds with mass greater 
than M > 1Q 5 5 M©, whose internal structure has a rea- 
sonable chance of being resolved, are plotted. (This 
cut-off is slightly lower (M > 1O 5 5 M ) than in TT09 
(M > 10 6 Mq) due to clouds being on average more 
massive). a c is the mass-averaged ID internal velocity 
dispersion of the clouds viewed at a 52° inclination an- 
gle, similar to our view of M33. 

Observed clo uds in M33 (iRosolowskv et al.ll2003D and 
the Milky Way (|Larsonlll981t ISolomon et al.lll987t ) show 
a velocity dispersion (line width) that increases as a 
power of their radii. The fit for these two observed popu- 
lations are also plotted in Figure GO Unlike in disk NoSF, 
where we reproduced the observed relation over an order 
of magnitude, neither of our disks have a cloud popula- 
tion that shows a strong correlation between these two 
quantities. This is likely due to the restricted range of 
radii that the clouds now obtain, covering only a factor 
of four in range, compared to in disk NoSF where they 
extended over a factor of 10 without the inclusion of star 
formation to limit growth. 

In part, this uniformity of our current cloud popula- 
tion radii is due to our initial conditions. We intention- 
ally model clouds in an environment typical to the Milky 
Way at the solar radius. The range in disk surface den- 



sity is only a factor of 10 as can be seen in Figure [2] and 
decreases over time. This means that any variation in 
cloud properties in our populations must come from their 
interactions and internal physics, not from variations in 
the galactic environment. The inclusion of energetic lo- 
calized feedback might be expected to ease this issue and 
create a wider range of cloud radii. 

The effect of diffuse heating in disk SF+PEheat 
slightly reduces the spread in cr c , but otherwise has no 
impact. This is in difference to the effect on the warm 
ISM, where the presence of diffuse heating significantly 
reduced the velocity dispersion, as seen in Figure [2] The 
denser, self-gravitating cloud gas, however, is less altered 
by this term. 

5.1.3. Simulation Time: Distribution of the Angular 
Momentum Vector 

The distribution of angles, 9, between the cloud's own 
angular momentum vector and the galactic rotation axis 
is shown in Figure [10] for both disks at four simulation 
times; 50, 100, 200 and 300 Myr. 9 for clouds in disk 
SFOnly are shown on the top panel while 9 for the clouds 
in Disk SF+PEheat are in the lower panel. Clouds with 
< 9 < 90 rotate about their center of mass in the same 
sense as the galaxy's angular momentum vector and are 
considered to be prograde rotators. Clouds with 90 < 
9 < 180 rotate in the opposite sense to the galaxy and 
are retrograde rotators. As the first clouds begin to form 
in the disk, we see an almost purely prograde population 
for both simulations. These clouds predominantly feel 
only the disk sheer as cloud-cloud interactions are still 
low in the partially fragmented disk. As the simulation 
continues, a retrograde population of clouds develops and 
by 300 Myr, 30% of the clouds in the disk without diffuse 
heating rotate in the opposite sense to the galaxy. Clouds 
present at these later times feel not only the sheer of 
the disk, but also the gravitational pull from neighboring 
clouds, including the effect of mergers. While forming in 
this more complex environment, clouds can develop a 
retrograde motion, or be later scattered by neighboring 
clouds to switch their sense of rotation. 

When diffuse heating is included in disk SF+PEheat, 
the fraction of clouds that develop a retrograde motion is 
significantly reduced, with only 12% rotating in a counter 
sense to the galaxy at 300 Myr. That this was true was 
indicated in Figure [5](e), where the specific angular mo- 
mentum profile contained only a small component with 
j z < 0.0. As mentioned in section IB.1.1[ this is the re- 
sult of the denser warm ISM, whose filamentary structure 
acts to encourage the clouds to remain prograde. 

The effect of t he ISM enviro nment on cloud properties 
was also noted bv lDobba (|2008h in their single fluid model 
of the ISM. In their disk simulation using an isothermal 
warm gas, no clouds arose with retrograde rotation, a 
fact they attribute to the absence of a clumpy medium 
required to cause collisional interactions between clouds. 
In our model, the ISM is fully multiphase, but the diffuse 
heating has increased the density and structure of the 
warm gas. We do not find a decrease in th e number of 
cloud mergers, but similarly to[Bobbs (200^), we find the 
warm ISM is playing a significant role in determining the 
cloud rotation. 

The result that the disk environment can impact the 
cloud properties is particular interesting, since that was 
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Fig. 9. — Cloud velocity dispersion vs size relation for clouds with mass M > 10 5,5 Mg at t = 200 Myr. Clouds in disk SFOnly are 
plotted in the left-hand pane and clouds in disk SF+PEheat are plotted on the right. a c is the line of sight velocity of the clouds in the 
galaxy viewed at 52° inclination angle, similar to our view of M33. The dashed blue line shows the result (slope and normalization) of the 
IRosolowskv et al.l | |2003|) st udy of massive (10fM Q < M c < 10 6 M ) GMCs in M33, with an exponent of 0.45 ± 0.02. The green dotted 
line shows the result of the Solomon et al. ( 1987) study of Galactic GMCs with an exponent of 0.5 ± 0.05. No obvious trend is seen in our 
results due to the limited range of radii of the clouds. 
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Fig. 10. — Distribution of the angle between the cloud angular momentum vector and the galactic rotation axis at different times over 
the course of the simulation. The distribution for clouds in disk SFOnly is shown in the top panel and in disk PEheat in the bottom panel. 
The shaded bars indicate retrograde motion. The presence of diffuse heating significantly decreases the number of retrograde clouds at late 
times. By 300 Myr, the number of retrograde clouds has reached 30% in the disk without diffuse heating and 12% in the disk with heating. 
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Fig. 11. — Distribution of cloud lifetimes in the simulations. 
Green circles show the clouds in disk SFOnly. Red squares depict 
clouds in disk SF+PEheat and black triangles show the results 
from disk NoSF. The majority of clouds live only a few Myr with 
only of order 0.5% reaching 100 Myr. This high end tail is probably 
due to the lack of energetic feedback. 

not a previous prediction from analytical models or ob- 
servations. In disk NoSF, the percentage of clouds that 
develop a retrograde motion was equal to disk SFOnly, 
supporting the claim that it is the difference in the warm 
ISM that controls t he development of retrog rade motion. 

Observations by iRosolowsky et al.l (|2003f ) find a 40% 
split between prograde and retrograde clouds in the sam- 
ple of GMCs in M33. Both our disks support the hypoth- 
esis that interactions between self-gravitating clouds are 
responsible for the distribution in angular rotation. 

5.2. GMC Properties with Cloud Age 

By tracking the clouds through the simulation over 
time, as described in section [2. 2 \ we are able to compare 
properties of clouds that have the same age, t', but exist 
at different simulation times. To avoid our results being 
affected by the initial fragmentation of the disk (during 
which cloud interactions are low) we only consider clouds 
born after 140 Myr in our analysis. 

The spread of ages in the clouds of our simulations are 
shown in Figure 1111 In the two simulations presented in 
this paper, the vast majority of clouds only live between 
''ago = — 20 Myr. This number drops by half after 
3 Myr and by a factor of ten by 20 Myr. Of the 7% that 
live longer than this, only 0.5% survive to 100 Myr. 

Without star formation, the decrease in cloud number 
is steady with age, as cloud destruction occurs purely 
as a function of the merger rate which is approximately 
constant. The addition of star formation increases the 
cloud mortality rate by a factor of 3-4 by a cloud age of 
10 Myr, having the greatest impact on the younger, less 
massive clouds (see below for details of the evolution in 
the cloud mass). It is likely these very young objects are 
forming close to other clouds who they almost immedi- 
ately merge with, or that their smaller size causes them 
to be disrupted quickly by star formation. 

Including the diffuse heating in disk SF+PEheat de- 
creases the star formation rate and allows clouds to live 
longer. It has the most effect on older clouds, whose 



higher mass makes them more resilient to mergers and 
hence their death is controlled primarily by their star 
formation. 

The lifetimes of GMCs is still a matter of intense de- 
bate. A current prevailing point of view, however, is 
that GMCs live 1-2 dyn amical times, putting their age in 
the r ange of 5-20 Myr ([McKee fc Ostrikerl [20071 iMurravl 
[20T1 . This agrees well with our clouds, a surprising re- 
sult since we include no form of energetic feedback which 
has been considered to be one of the main factors in cloud 
destruction. 

5.2.1. Cloud Time: Cloud Property Distributions 

The properties of the clouds as a function of their age 
are shown in Figures [T^] (disk SFOnly) and Q2] (disk 
SF+PEheat) for the same distributions plotted in Fig- 
ure [7] and [H cloud mass, radii, surface density, virial 
parameter, vertical component of the angular momen- 
tum and their vertical distribution. The solid black line 
shows recently formed clouds that are less than 1 Myr 
old. The red dashed line shows the distributions for 
clouds between 9-10 Myr old, the estimated observed age 
for a GMC. The blue dot-dashed line is for cloud ages 49- 
50 Myr and the green dotted line for clouds 99-100 Myr 
old. These two last lines show the properties of the older 
clouds in the simulation although, as shown in Figure [TT1 
significantly fewer clouds live to these ages, with only 44 
clouds forming the distribution of the 99-100 Myr profiles 
in disk SFOnly and 92 clouds in disk SF+PEheat. 

In both the disks, we see that older clouds are on aver- 
age more massive, have larger radii, higher surface den- 
sity and are slightly more gravitationally bound. These 
trends are more pronounced than in disk NoSF since the 
gas depletion from star formation destroys a large num- 
ber of younger clouds, as seen in Figure [IT] Clouds with 
masses M < 10 5 are particularly prone to early mortal- 
ity through star formation and mergers. The clouds that 
survive their early years are more massive, accumulating 
gas through accretion and mergers at a faster rate than 
their star formation can deplete it. This effect is most 
marked in disk SFOnly. In disk SF+PEheat when diffuse 
heating is present, less clouds with masses < 10 5 M are 
born due to the additional heating increasing the Jeans 
length, raising the size of an object formed through grav- 
itational collapse. This reduces the evolution of the low 
mass distribution at in the first 10 Myr of the cloud's life. 

In contrast to this, the distributions of vertical angular 
momentum, Figure [l2](e) and Figure [T3](e), and vertical 
position above the disk, Figure IT2"1(f ) and Figure [131(f). 
are largely independent of cloud age. Without diffuse 
heating, newly formed clouds in disk SFOnly have 
a very low vertical angular momentum, causing the 
distribution to be sharply peaked. Cloud interactions 
later broaden this profile slightly and reduce the peak, 
but the form remains largely unchanged. When diffuse 
heating is present in disk SF+PEheat, the clouds are 
born with a wider range of positive (prograde) angular 
momentum values that also show little change over time. 
The filamentary structure of the ISM has produced a 
more dominantly prograde cloud population as seen 
in Figure 1101 and the higher values of j z means the 
probability of later producing a retrograde cloud in a 
cloud collision is low. 
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Fig. 12. — Normalized distributions of the GMC properties described in Figure [7] for disk SFOnly, but now showing results for different 
cloud ages: 0-1 Myr (solid lines), 9-10 Myr (dashed lines), 49-50 Myr (dot-dashed lines) and 99-100 Myr (dotted lines). All clouds in this 
analysis were born after 140 Myr of disk evolution, i.e. in the fully fragmented phase. Top left, (a): cloud mass, M c . Middle left, (b): cloud 
radius, Rc,A- Bottom left, (c): mass surface density, E c . Top right, (d): virial parameter, a v ir- Middle right, (e): vertical component of 
the specific angular momentum, j z . Bottom right, (f): cloud center-of-mass vertical positions, z. 
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Fig. 13. — Normalized distributions of the GMC properties, as described in Figure [51 for Disk SF+PEheat, for cloud ages: 0-1 Myr (solid 
lines), 9-10 Myr (dashed lines), 49-50 Myr (dot-dashed lines) and 99-100 Myr (dotted lines) for clouds born after 140 Myr of disk evolution, 
i.e. in the fully fragmented phase. Top left, (a): cloud mass, M c . Middle left, (b): cloud radius, R C ,A = (Ac/tt) 1 / 2 . Bottom left, (c): 
mass surface density, S c . Top right, (d): virial parameter, o v i r . Middle right, (e): vertical component of the specific angular momentum, 
j z . Bottom right, (f): cloud center-of-mass vertical positions, z. 
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5.2.2. Cloud Time: Distribution of the Angular Momentum 

Vector 

The robustness of the angular momentum with respect 
to cloud age is shown in the distribution of angles, 8, be- 
tween the cloud angular momentum vector and the galac- 
tic rotation axis plotted in Figure Q31 Clouds at four dif- 
ferent ages are shown; newly born clouds with t' < 1 Myr, 
then clouds with ages t' = 9—10 Myr and 49—50 Myr and 
finally our oldest clouds with t' = 99 - 100 Myr. The top 
panel shows the distribution of for disk SFOnly while 
the lower panel shows the distribution for clouds in disk 
SF+PEheat. As with Figure [TOl the shaded bars show 
angles that equate to a retrograde rotation. 

The percentage of retrograde rotating clouds remains 
about 25% for all age ranges in disk SFOnly and about 
10% in disk SF+PEheat. Since all these clouds are born 
after 140 Myr in the simulation, they feel the forces from 
the fully fragmented disk, in addition to the sheer. The 
gravitational interactions from neighboring clouds on the 
newly forming bodies cause the retrograde population to 
form. As we can see from the first panel of Figure [TU1 
clouds forming in pristine gas are always prograde. 

As we have seen previously, the effect of diffuse heating 
is to reduce the fraction of retrograde clouds, due to the 
filamentary structure of the warm ISM having a strong 
impact on the newly forming clouds. 

6. STAR FORMATION 

The star formation history for the disk is plotted in 
Figure 1151 As with the cloud formation history in Fig- 
ure [5J the star formation initially increases steeply as the 
disk fragments. It reaches a peak value around 125 Myr 
and then declines steadily as gas is consumed in the disk. 

At 200 Myr, the star formation rate is ISMQyr" 1 for 
disk SFOnly and 17 Meyr" 1 for disk SF+PEheat. By the 
end of the simulation, this has dropped to 3M Q yr -1 and 
lOMQyr^ 1 , respectively, although the lower numbers are 
purely a factor of the gas depletion. The Milky Way 
is estimate d to have a star formation rate of order 1 — 
^M^ yr" 1 ([Murray fc Rahman|[2010HWilliams fe McKed 
1997). The fact we are higher than this though, is not 
surprising due to our lack of localized feedback. 

The addition of diffuse heating in disk SF+PEheat ini- 
tially reduces the fraction of dense cloud gas in the disk, 
as we saw in Section [4] causing the SFR to be lower 
over the first ~ 175 Myr than when heating was not in- 
cluded. This causes the gas to be depleted in the disk 
at a slower rate, resulting in more gas being available at 
later times. This can be seen in Figure [161 which shows 
the distribution of the fraction of mass in stars for the 
clouds present at three different simulation times. Early 
in the simulation, clouds are gas dominated with a low 
(< 0.1) stellar fraction. By 300 Myr, a large fraction of 
the gas has converted into stars, making the majority of 
the clouds strongly stellar dominated. This is especially 
true for clouds in disk SFOnly, where 3/4 of the clouds 
have a stellar fraction > 0.9 by the end of the simulation. 
This inevitably causes the star formation to decrease and 
after 175 Myr, the star formation rate in disk SFOnly has 
dropped below the rate in disk SF+PEheat because the 
gas abundance has become too low. 



The Kennicutt-Schmidt relation is an empirical mea- 
surement of the relationship between the surface SFR, 
S s f r , and the surface gas density, £ ga s, in a galaxy. It 
describes how efficiently a galaxy is converting gas into 
stars and takes the form shown in Equation [3] The 
value of the expo nent, q s fr, has b e en measured by sev- 
eral groups J^ lKennicutti[l998l: I Wong fc Blitl 120021: 
IBigiel et all I2008D for both the Milky Way and other 
galaxies and found it to vary between a s f r ~ 1 — 3. Most 
recent work (|Bigiel et al.ll2008t IWong fc Blitl2002ft sug- 
gests that the correlation is truly with molecular gas, 
Eh 2J rather than total gas, S gas = Ehi + £h~ 2 , f° r which 
Qsfr,H 2 = 1-0 + 0.2. When averaged with the atomic com- 
ponent, the exponent increases before the relation breaks 
down around S gas ~ IOMqpc" 2 in the outer regions of 
the disk where the gas becomes saturated with HI. 

Figure [T7] shows this relation plotted for our disks at 
t = 200 Myr. The diagonal lines mark constant star for- 
mation efficiency and the black crosses show the obser- 
yations from The HI Nearby Galaxy Survey (THINGS) 
(|Bigiel et al.l [2008). For these observations, the fall-off 
around 10 M Q pc~ 2 is due to the dominance of atomic gas, 
with the points at higher densities with a s f r = 1.0 ± 0.2 
being purely from molecular i?2- The results for disk 
SFOnly are shown by green filled symbols while disk 
SF+PEheat is marked in red open symbols. Square sym- 
bols show the result of averaging over an area of 750 pc 
across around each cloud, the equivalent spatial resolu- 
tion to the THINGS survey. The circles show the quanti- 
ties averaged over the individual clouds themselves, with 
area of order 15 pc across. 

For equivalent spatial resolution (squares) we see that 
our star formation rate is higher than the observations 
by a factor of 10 in both our simulations. This agrees 
with what we saw in Figure [T5] and is likely due to not 
having a source of localized feedback. The gradient here 
is a s f r 1.77 for the disk without diffuse heating and 1.81 
for when heating is included. This is steeper t h an th e 
fit for the molecular gas found by IBigiel et~aT1 (|2008h . 
but in closer agreement with the gradient of the atomic 
gas. Given the size of the region we are averaging over, 
compared to our cloud size, we would expect our result 
to be atomic gas dominated. 

When averaged over individual clouds, the gas surface 
density is higher by approximately a factor of 100, com- 
pared with averaging over a larger area. In this volume, 
we can assume that the gas is at least 50% molecular, 
depending on the mass of the assumed atomic envelope. 
The fitted gradient for both these cloud p opulations is 
1.27, slightly steeper than the result from ([Bigiel et al.l 
2008), likely due to this mix of atomic and molecular gas. 

Note, that neither sets of points gets a gradient of 1.5, 
as would be expected if our results were simply a product 
of having a constant star formation efficiency term, as 
given in Equation O 

Despite having a lower SFR over the first 2/3 rds of the 
simulation, we do not see any great differences between 
the cloud populations with and without diffuse heating. 
From Figure [15l we can see that at our time of analysis, 
t = 200 Myr, the SFRs in both simulations are approxi- 
mately constant. However, the same plot at t = 150 Myr 
shows no greater disparity between them, due their dif- 
ference in SFRs being accompanied by similar differences 
in gas surface density, as was seen in Section [5.11 There 



21 



io J 





| 


1 1 1 1 




| 


1 1 1 1 




I 


1 1 1 


1 1 


| 


1 1 1 


1 


r 
r 






- 1 Myr " " 






9- 10 Myr " 

:: 

T 






. siiP - : ■ 


49 -50 Myr " 






99- 100 Myr" 




1 






1 






i 






1 








45 


90 135 180 
6 [deg] 


45 


90 13 
9 [deg] 


^ 180 


45 


90 135 180 
6 [deg] 


45 


90 1 
G [deg] 





45 90 135 180 45 90 135 180 45 90 135 180 45 90 135 180 

[deg] [deg] [deg] [deg] 



Fig. 14. — Distribution of the angle between the cloud angular momentum vector and the galactic rotation axis for clouds of different ages, 
born after 140 Myr, in disk SFOnly (top) and disk SF+PEheat (bottom). The shaded bars indicate retrograde motion. The percentage of 
clouds rotating retrograde in the top row is 28%, 26%, 18%, 25% respectively for each panel and 9%, 10%, 8%, 4% for the bottom row. 



is slightly less scatter in the simulation that includes the 
heating term, especially when averaged over 750 2 pc 2 , 
due to the less fragmented ISM distribution. 

Focusing on the individual clouds, the specific star 
formation rate (SSFR; the star formation rate per unit 
mass) can be plotted as a fu nction of th e cloud 's gas 
mass. In their recent paper, lLada et all (|2010D mea- 
sure the SFR in eleven GMCs in the Milky Way. They 
find a linear dependence with the cloud's mass, specif- 
ically SSFR= 4.6 ± 2.6 x lO^yr" 1 . The relation for 
our clouds is plotted in Figure [15] for the SSFR aver- 
aged over one dynamical time. The left-hand plot shows 
the results for disk SFOnly and the right-hand plot for 
disk SF+PEheat. Black circles show binned data for the 
complete cloud sample while the red squares and blue tri- 
angles only include clouds that have undergone a merger 
in the last dynamical time, with the latter considering 
only major mergers with a cloud mass ratio less than 
2.0. 

The upper dotted line shows the value from lLada et al.l 
(2010). With the exception of the most massive clouds, 
our s tar formation rate is lower by a factor of 3-4. How- 
ever, |Lad^eOI] (|2010[) calculate the cloud mass from gas 
above a threshold surface density of E gas « 116 M pc -2 , 
which they estimate is equivalent to a volume density 
of tih 2 ~ 10 4 cm~ 3 . While our clouds larger than 
M « 10 5 Mq achieve this surface density, their volume 
density is much lower, with the majority of clouds having 
ave rage densities betw een 10 2 — 10 3 cm~ 3 . If we reduce 
the lLada et al.l (|2010h SSFR value by a factor of 10 to 



allow for this, we get the lower dotted line. Our cloud 
SSFR largely lie above this lower line, suggesting that 
the SFR per cloud is too high, given the resolution, by 
a factor of 10 in agreement with Figure [IT] for the disk 
averages over the same spatial scale as the observations. 

For clouds with masses M < 10 5 5 M Q , the SSFR is 
approximately constant. At first sight, t his seems to 
be in agreement with the lLada et al.l (|2010l ) result, how- 
ever this population of clouds are at our resolution limit 
where their internal structure cannot be resolved well. In 
our cloud population, density scales with mass so larger 
objects have shorter dynamical times, increasing their 
SSFR. In this region where our results are not resolution 
limited, M > 1O 5 M , we therefore find that the SSFR 
is proportional to the cloud mass. 

This result remains when we introduce diffuse heating 
in the right-hand plot of Figure[T8] The maximum SSFR 
however is less, in agreement with our previous findings, 
and there is evidence of a down turn at high masses. This 
small population of very massive (> 1O 6 5 M ) clouds 
are extended structures from recent mergers, resulting in 
them having a lower density. 

As we saw in Figure [5] mergers between clouds are a 
common occurrence and likel y to have a significant im- 
pact on the clouds evolution. iTanl (|2000h suggested that 
such collisions could trigger star formation, providing 
a way of connecting local-scale motions with the glob- 
ally observed Kennicutt-Schmidt relation. In Figure [T8l 
however, we find that the presence of a recent merger de- 
creases the clouds SSFR. This effect is even greater if the 
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Fig. 15. — The star formation rate over the course of the simula- 
tion for disk SFOnly (green solid line) and SF+PEheat (red dashed 
line). Like the cloud population, the SFR peaks shortly after the 
disk has become fully fragmented. As gas is converted into stars 
and depleted from the clouds, it begins to drop. This decrease in 
SFR is significantly slower when diffuse heating is included, drop- 
ping to a value of 10 Mq yr — 1 rather than 3Mq yr — 1 when diffuse 
heating is not included. Around 200 Myr, before gas depletion has 
removed most of the gas in the clouds with no diffuse heating, the 
SFR is greater than the estimated Milky Way value by a factor of 
10. This is probably due to the lack of localized feedback. 



merger was a major one with a cloud whose mass was at 
least 50% of the object it merged with. The impact of the 
mergers depends on mass, with larger clouds being un- 
affected unless experiencing a major merger. This is not 
hard to understand, since a more massive cloud will be 
undergoing predominantly collisions with much smaller 
objects which will cause less perturbations to its inter- 
nal structure. The reason for the drop in SSFR is due 
to a corresponding drop in the cloud density. However, 
such a result is liable to be resolution dependent, due to 
the difficulties in resolving the changes t o the clouds' in- 
ternal structure across only a few cells. iFederrath et al.l 
(2010) finds that 30 cells are required to resolve a vortex 
within a cloud, while local box simulations with resolu- 
tions on the sub-parsec scale suggest that low density 
(uh = 3 cm -3 ) colliding flows can trigger local star for- 
matio n by initiating gravitational collapse (jHeitsch et al.1 
2008). If this applies to higher density collisions, then it 
is probable that cloud collisions could help, rather than 
hinder, star formation. 

7. CONCLUSIONS 

We investigated the formation and evolution of the 
GMC population formed in two isolated, Milky Way- 
type galaxies at a resolution < 10 pc. Both our simu- 
lated disks included star formation and radiative cooling, 
with our second model also including a diffuse heating 
term, representative of the photoelectric heating from 
dust grains. We did not include any form of localized 
energetic feedback (e.g. supernovae). 

Both the disks fragmented through gravitational insta- 
bilities to form a population of clouds, which we identi- 



fied as GMCs when their densities reached the threshold 
value of 7t,h,c > 100 cm~ 3 . These were tracked over time 
to provide statistics both as a function of simulation time 
and of cloud age. 

The number of clouds formed above M c > 10 5 M Q 
was found to agree well with the observed Milky Way 
GMC population. The properties of the clouds were also 
comparable to observations, including the distributions 
of cloud mass, size, mass surface density, virial param- 
eter, angular momentum, vertical height above the disk 
and the distribution of angles of angular momentum with 
respect to the galactic rotation axis. 

Cloud ages were found to lie largely between 0—20 Myr, 
in good agreement with current estimates. It is no- 
table that this is without a source of localized feedback 
which has previously been expected to provide a dom- 
inant mechanism for cloud destruction. Many of our 
clouds die in the first 3 Myrs, reducing our population 
by 50%. This cloud infant mortality is due to mergers 
with nearby forming clouds and star formation which can 
destroy low mass (M < 1O 5 M ) objects. 

The inclusion of diffuse heating raised the pressure of 
the warm and cold ISM to suppress the fragmentation 
of the disk, producing an initially smaller population of 
clouds embedded in a more massive and structured warm 
ISM. The denser ISM retains a filamentary structure af- 
ter fragmentation that maintains a lower velocity disper- 
sion than when heating was absent. This environment 
has two major impacts on the cloud properties: 

1. The filamentary warm ISM produces a predomi- 
nantly prograde rotating population of clouds even 
at late times. Without diffuse heating, the clouds 
become 1/3 rd retrograde both with and without 
star formation. This is evidence that the environ- 
ment of the cloud plays a dominant role in its evo- 
lution. 

2. The second effect is that the lower mass of cloud 
material in the disk reduces the star formation rate 
in the first 175 Myr of the simulation. Past this 
time, the star formation rate remained approxi- 
mately constant, while in the absence of heating, 
gas depletion in the clouds causes the production 
of stars to drop by a factor of 10. 

Cloud mergers and interactions were a frequent occur- 
rence in both disks, occurring at a rate of ~ 0.25 of 
an orbital period. This is only slightly higher than the 
merger rate recorded in TT09 of the same simulation 
without star formation. The effect of mergers appears 
to be to reduce the cloud density, thereby reducing the 
star formation rate, but we note this effect is likely to be 
dependent on resolution. 

The star formation rate in the disks was roughly a 
factor of 10 too high in both the disk without diffuse 
heating and when it was included. There are several 
possibilities why this could be the case including the ab- 
sence of localized feedback and processes such as added 
support from magnetic fields. Previous work that ex- 
plores the impact of magnetic fields on the GMC popu- 
lation suggests that the magnetic pressure can suppress 
the format ion of the popu l ation. However, research per- 
formed bv lDobbs fc Price! (|2008t ). suggests that the gas 
to magnetic pressure must be (3 < 0.1 for the magnetic 



23 



t [Myr] 
100 
200 
300 



g 0.4 ■ 

z° 



0.4 0.6 0.8 

M /M , 



g 0.4 

z" 



t [Myr] 
100 
200 
300 



"1-1 



J 



0.4 0.6 0.8 

M /M , 



Fig. 16. — Fraction of mass in stars for clouds at simulation times 100 Myr (solid blue line), 200 Myr (dashed black line) and 300 Myr 
(red dot-dashed line). Left plot shows results for disk SFOnly while the right-hand plot shows the clouds in disk SF+PEheat. 



fields to suppress fragmentation in cold gas. Alone, the 
magnetic Parker instability is not thought to be able 
to produce a realistic population of GMCs, although it 
might play a role in s eeding the gravitational collapse 
(|Kim fc Ostriker|[200l) . At the other end of the scale, 
localized feedback from sources such as supernovae or ra- 
diation driven winds can destroy clouds once they have 
formed stars. The exploration of both these mechanisms 
will be explored in future studies. 

In addition to the physical processes not yet included, 
a second source of error has to be the limit of our resolu- 
tion. On average, our clouds contain 76 cells and, for an 
average radius of 16 pc, this equates to roughly 4 cells 
in each dimension, although it is worth noting that our 
average radius in the plane of the disk is larger, at 34 pc. 
The impact this has on our calculations for the cloud 
properties was investigated in TT09. Figure 11 in that 
paper shows the effect of reducing our limiting resolution 
to 15.6 pc and 31.2 pc, plotting the cloud properties 
shown in Figures [7] and [8] for the populations in these 
runs. We found that clouds with masses M > 10 6 M Q 
were converged at all resolutions, but the peak mass, 



radius and virial parameter decreased with increasing 
refinement. Of particular concern was the sensitivity 
of the rotationa l dyna mics of the clouds to resolution. 
iFederrath et~aTl (|2010f l find that 30 cells are needed to 
accurately resolve a vortex, which is almost a factor of 
ten more than the resolution we are currently able to 
achieve within our clouds. To test the impact of this, we 
plotted the variation of the clouds' angular momentum 
with respect to the disk's rotation in Figure [14] for 
clouds with masses greater than M > 10 6 M Q , but saw 
no change in the distribution of 6. This is a promising 
indication that our result will remain true for a more 
highly resolved cloud structure, but we note that we are 
not able to test this directly at the present time. 

The author would like to thank Jonathan Tan, Ralph 
Pudritz and James Wadsley for helpful discussions. 
EJT also acknowledges the University of Florida High- 
Performance Computing Center for providing computa- 
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Fig. 18. — Specific star formation rate averaged over one dynamical time (cloud free-fall time) vs. cloud gas mass for individual clouds 
in disk SFOnly (left plot) and disk SF+PEheat (right). Black circles show the distribution for all clouds, red squares include only clouds 
that have undergone a merger in the last dynamical time and blue diamonds include only the clouds for whom that merger was with a 
cloud greater than 50% of its mass. 



